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Abstract 



We investigate the chiral properties of quenched domain-wall QCD 
(DWQCD) at the lattice spacings ~ 1 and 2 GeV for both plaquette 
and renormalization-group (RG) improved gauge actions. In the case of the 
plaquette action we find that the quark mass defined through the axial Ward- 
Takahashi identity remains non-vanishing in the DWQCD chiral limit that 
the bare quark mass mj ^ and the length of the fifth dimension A'^, ^ oo, 
indicating that chiral symmetry is not realized with quenched DWQCD up to 
a^^ ~ 2 GeV. The behavior is much improved for the RG-improved gauge ac- 
tion: while a non-vanishing quark mass remains in the chiral limit at a^^ ~ 1 
GeV, the result at ~ 2 GeV is consistent with an exponentially vanishing 
quark mass in the DWQCD chiral limit, indicating the realization of exact 
chiral symmetry. An interpretation and implications are briefly discussed. 
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I. INTRODUCTION 



A recent exciting development in lattice field theory is the realization of exact chiral-like 
symmetry at finite lattice spacing without fermion doubling. It started with the domain- 
wall and the overlap 0] formalisms, which were originally proposed to formulate chiral 
gauge theories. They were soon applied to vector gauge theories such as QCD and 
progress culminated in the recent discovery |^ that the essential feature of these formalisms 
is the existence of chiral-like symmetry which follows from the Ginsparg- Wilson relation . 
An explicit form of the fermion action satisfying this relation is already known 0]. These 
developments lead us to expect that QCD with exact chiral symmetry can be numerically 
simulated on a lattice. 

The domain-wall fermion formalism is a five dimensional extension of the Wilson fermion 
action with a negative mass —M. The theoretical basis is that the effective four- dimensional 
theory obtained by integrating out the heavy unphysical modes satisfies the Ginsparg- Wilson 
relation [^,|Tl|]. On the practical side, the knowledge accumulated on the Wilson fermion 



action over the past twenty years enables an efficient implementation of this system on 
computers. 

There are, however, two subtleties that need to be clarified for numerical applications of 
the domain wall formalism to QCD (DWQCD). First, the length of the fifth dimension, A^^, 
cannot be set to infinity: one has to study the finite A^^ effect. Second, one has to tune the 
domain- wall height (i.e., fifth dimensional mass) M to an appropriate value in order to keep 
the massless mode. 

In the free theory one massless mode exists in the range < M < 2 0]. The two 
chiral modes, with opposite chiralities, are exponentially bound to the opposite boundaries 
in the fifth dimension. The situation remains unchanged in perturbation theory. An explicit 



calculation at one-loop order shows [12,13| the existence of the massless mode guaranteed 



for the range < M' < 2 where M' receives a shift by the one-loop quantum correction. It 
is also easy to demonstrate that finite A^^ effects are suppressed exponentially in A^. 

Numerical tests of chiral properties and physical applications of (quenched) DWQCD 
were first discussed in Ref. Through measurements of pion mass and the K^—K^ mixing 



matrix element at the gauge couphng 13 = Q/g^ 6.0 (a^^ ~ 2GeV), these quantities were 
shown to vanish in the chiral limit, as required by chiral symmetry, for a rather small fifth 
dimensional length of A^s ~ 10. Good chiral behavior has also been observed for A^s ~ 10 



for a gauge field background with non-trivial topology [15]. Even the Atiyah-Singer index 



theorem has been shown to be approximately satisfied on a lattice for a non-zero quark mass 



15,16 



These successful results have promoted further applications. The calculations include the 
strange quark mass using nonperturbative renormalization factors |r^l9|, the pion decay 
constant ||20|, the hadron mass spectrum |2l|j2^ and QCD thermodynamics with dynamical 
quarks p3|-p5[ 



Earlier studies with DWQCD were limited to lattice spacings of around ^ 2 GeV. 
It is obviously necessary to investigate the scaling behavior and carry out the continuum 
extrapolation. Reducing the lattice spacing from ~ 2 GeV, however, is computationally 
difficult, since the formalism requires Ns times more CPU time than usual simulations. 
Simulations at coarser lattices are easier, but we must clarify how large one can increase the 
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lattice spacing while keeping the massless mode. 

The chiral properties of quenched DWQCD have been examined in a number of recent 
reports p2| , |26| -p^ employing large fifth dimensional lengths up to Ng ~ 50. The authors 



of Refs. [p2| , p6| , who made a study with the plaquette gauge action at (3 = 5.7 and 5.85 
(a^^ ^1 — 1.5 GeV), claim that pion mass at such strong couplings does not vanish in the 
chiral limit after the Ng ^ oo extrapolation, indicating that exact chiral symmetry is not 
realized. Our own study at /5 = 5.65 [E3] yields the same conclusion. We also find that a 



renormalizat ion-group (RG) improved gluon action ||29| does not improve the situation at 
a-^ ~ 1 GeV 



In this paper we present a more systematic study of the chiral property of DWQCD 
covering both strong and 'weak' coupling regions, corresponding to lattice spacings a^^ ~ 1 
and 2 GeV with both plaquette and RG improved actions. We emphasize that the chiral 
properties are best studied with the quark mass extracted by the Ward-Takahashi (WT) 
identity for the axial vector current. The difference of this mass and the bare quark mass 
defines what we call 'anomalous quark mass'. We discuss that this quantity suffers from 
much less systematic uncertainties than the pion mass from both theoretical and numerical 
points of view. We study the behaviors of the anomalous quark mass as a function of A^^ and 
examine whether it vanishes exponentially in the Ng ^ oo limit for various coupling regimes, 
as required by the self-consistency of the proper chiral theory. This study is also made with 
varying M and a four dimensional spatial extent at some fixed coupling. We carry out this 
study with a parallel set of simulations employing both plaquette and RG-improved gauge 
actions to see the effect of improvement of the gauge field dynamics on the chiral property. 
In connection to our work we refer to Ref. which attempts to extract the WT identity 
quark mass from the pseudoscalar susceptibility. 

This paper is organized as follows. In section |I| we define the action and the anomalous 
quark mass. Numerical simulations and run parameters are described in section |T|. In 
section ^ we investigate the chiral property for the plaquette action, and in section |V] for 
the RG-improved gauge action. Our provisional interpretation is given in section for 
the results we found from the simulation. We close the paper with a brief summary and 
comments on the application of DWQCD (section |VI1| ). Appendix A presents an analysis 
with pion mass, and discuss whether they support the conclusions based on the anomalous 
quark mass. Numerical data for the anomalous quark mass and pion mass are collected in 
Appendix B. 



II. ACTION AND AXIAL WARD-TAKAHASHI IDENTITY 



We adopt Shamir's domain- wall fermion action @J§], except that the Wilson term and 
the domain wall height M have minus signs: 



Sf = - J2 i^{XiS)Dd-u,f{x,S]y,s')i){y,s') +^mfq{x)q{x) , 

x,s,y,s' X 

Ddwf{x, s] y, s') = D^{x, y)Ss,s' + D^{s, + (M - 5)6cc,ySs,s' , 

D\X, = E ^ [(1 - l,)UxJx-,f.,y + (1 + 7,)UlJx-i., 



(II.l) 

(11.2) 
(11.3) 
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PlS2,s' {s = 1) 

D%s,s') = { PlSs+i,s' +PrSs-i,s' {l<s<N,) , (II.4) 

Pr6n,~i,s' {s = Ns) 

where x, y are four- dimensional space-time coordinates, and s, s' are fifth-dimensional or 
"flavor" index, bounded as 1 < s, s' < Ng with the free boundary condition at both ends 
(we assume Ng to be even); Pr/l is the projection matrix Pr/l = (1 ± 75)/2, and m/ is 
physical quark mass. 

In DWQCD the zero mode of a domain-wall fermion is extracted by the "physical" quark 
field defined on the edges of the fifth dimensional space, 

q{x) = PLtpix, 1) + Pr.^{x, Ng), 

q{x) = ^{x, Ns)Pl + i^ix, 1)Pr. (II.5) 

The QCD operators are constructed from these quark fields. 

For the gauge part of the action we employ the following form in 4 dimensions: 



5'giuon = \{co Tr[/p, + ci Yl Trf/rtg + C2 Trf/chr + C3 Y Trf/p,g 

3 I plaquette rectangle chair parallelogram 



(11.6) 



where the first term represents the standard plaquette action, and the remaining terms are 
six-link loops formed by a 1 x 2 rectangle, a bent 1x2 rectangle (chair) and a 3-dimensional 
parallelogram. The coefficients Cq, ■ ■ ■ , C3 satisfy the normalization condition 

Co + 8ci + I6C2 + 8C3 = 1. (II. 7) 

The standard Wilson plaquette action is given by Cq = 1, Ci = C2 = C3 = . The RG- 
improved action of Iwasaki is defined by setting the parameters to Cq = 3.648, Ci = 
—0.331, C2 = C3 = 0: with this choice of parameters the action is expected to realize smooth 
gauge field fluctuations approximating those in the continuum limit better than with the 
unimproved plaquette action. 

The chiral transformation of domain wall fermion is defined as 0, 

dip{x, s) = iQ{s)e''T''-^{x, s), (II.8) 
5i^{x, s) = -ii^{x, s)Q{s)e''T'', (II.9) 

where Q{s) = sign(A^s — 2s + 1), T° is a generator of flavor symmetry group SU(A^/) with 
a transformation parameter. The corresponding WT identity becomes 

Y (V,A;:(a;)0) = 2mf {P\x)0) + 2 {jt,{x)0) + i {6^0) , (11.10) 
where the axial- vector current ^^(a;) and the pseudoscalar density P'^{x) are defined by 
K^x) = Y Q(s)^(^(^, s)T''{l - l,)U,{x)i,ix + s) 

s ^ 

-i^{x + /i, s){l + ^,)Ul{x)Ty{x, s)) , (11.11) 
P\x) = q{x)T'^^^q{x), (11.12) 
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and J^q{x) representing explicit breaking of chiral symmetry takes the form, 

W = -^i^, j)T'^PlHx, y + 1) y + 1)T'^Pr^{x, y). (11.13) 

For a smooth gauge field background the anomalous contribution (^J^q{x)0'^{y)^ is 
bounded by an exponentially small value with the argument —Ng and vanishes in the limit 
A^s oo 0. In this paper we adopt the operator O = P^(0,0), and measure the chiral 
symmetry breaking effect by 



"^« = .1-E,(P.(«.X)P'(0.0)>' '"-l^) 

which we call 'anomalous quark mass'. We expect this quantity to be dominated by short- 
range fluctuations since it is determined by the coupling strength of pion to the operators 
and P*. Therefore, finite-size effects are likely to be small. For the same reason, there 
should be no quenched chiral singularities in the relation between m^q and the bare quark 
mass rrif. 

The axial WT identity is written in terms of the anomalous quark mass as. 



E,V,^;(i,x)P>(y), 

The left-hand side of this equation, proportional to the pion mass squared, vanishes at 
mj + = 0. 



III. RUN PARAMETERS AND MEASUREMENTS 

We explore the chiral properties of DWQCD at both a^^ ^ 1 and 2 GeV using the 
plaquette and the RG-improved actions. For the plaquette action, we choose the couplings 
P = 5.65 and 6.0, which correspond to = 1.00 and 2.00 GeV, respectively, if one 
determines the scale from the string tension a |30| assuming cr = (440 MeV)^. For the 



RG-improved action, we take P = 2.2 and 2.6 for which a = 0.97 and 1.94 GeV also from 



the string tension reported in Refs. pT|-p5 



We generate quenched gauge configurations on four-dimensional A^"^ x A^^ lattices. One 
sweep of the gauge update contains one pseudo-heatbath and four overrelaxation steps. For 
runs at ~ 1 GeV we adopt the lattice size of 12^ x 24 and 16^ x 24 to study finite- 
size effects. At a^^ ~ 2 GeV only a 16^ x 32 lattice is employed. After a thermalization 
of 500 (2000) sweeps hadron propagators are calculated at every 100th (200th) sweeps for 
~ 1 GeV (2 GeV). 

The domain-wall quark propagator is calculated on A*"^ x Nr x Ng lattices, where the 
gauge configuration on each fifth dimensional coordinate s is identical and is fixed to the 
Coulomb gauge. The size of the fifth dimension is varied from A^^ = 10 to 50 at ~ 1 GeV, 
and from A^^ = 4 to 24 — 30 at ~ 2 GeV, depending on the gauge action. 
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To allow a chiral extrapolation, we employ fermion masses of ruf = 0.03,0.05,0.1 (at 
a-^ ~ 1 GeV) and mj = 0.02,0.04,0.06 (at a'^ ~ 2 GeV). The domain-wall height de- 
pendence is examined at ~ 1 GeV with the choice of M = 1.3, 1.7,2.1, and 2.5. At 
a'^ ~ 2 GeV we choose M = 1.8. 

The propagators are calculated with the conjugate gradient algorithm with an even-odd 
pre-conditioning. We place the source at s = 1 and A^^ in the fifth coordinate, so that 
the propagator from the physical field q to the domain-wall field ilJx,s for arbitrary (x, s) 
is obtained. Both local and exponentially smeared sources are employed in the spatial 
directions, and meson masses are measured for all possible combinations of quark masses 
while only degenerate combinations are evaluated for anomalous quark masses. 

Our simulation parameters and the number of configurations are given in Table |. We 
note that gauge configurations are generated independently for each choice of Ng and M, 
i.e., there are no correlations between the data generated with a different domain- wall height 
and a fifth-dimensional size. 

In Figs. and H we show typical data for the ratio of two-point functions defined in ([1I.14|) 
as a function of the temporal distance t. We obtain the anomalous quark mass m^g at each 
mf,M and Ng, by fitting the plateau with a constant, the fitting range being determined 
by the inspection of plots for the ratio and those for the effective pion mass. We choose 
8 < t < 16 as the fitting range for all simulations at /5 = 5.65 (plaquette) and (3 = 2.2 (RG). 
We use 10 < t < 22 at /5 = 6.0 for the plaquette action and /5 = 2.6 for the RG action. The 
numerical values for m^q and pion mass are given in Appendix B. 

IV. PLAQUETTE GAUGE ACTION 

We discuss the anomalous quark mass for the plaquette gauge action. The bare quark 
mass {rrif) dependence of the anomalous quark mass m^g is illustrated in Fig. ^for /3 = 5.65 
and 6.0. We observe only a mild dependence of m^g on mj. The anomalous quark mass in 
the nominal chiral limit m^g{mf = 0,Ns,M) is obtained by a linear extrapolation in nif. 
The extrapolated values are collected in Appendix B. 

If chiral symmetry is realized in the Ng ^ oo limit of the DWQCD system, m^g{mf = 
0, Ns, M) should vanish exponentially in Ng. We examine this point first at a strong coupling 
(3 = 5.65 corresponding to ~ 1 GeV. The anomalous quark mass in the nominal chiral 
limit, m^q{mf = 0,Ns,M), is plotted as a function of A^^ in Fig. ^ for three values of M 
(M = 1.3, 1.7, 2.1). A comparison of the filled symbols taken on a 16^ x 24 x lattice with 
the open ones on a spatially smaller 12^ x 24 x A^^ lattice shows no significant finite-size 
effect between the two spatial sizes corresponding to A'^-a ~ 2 and 3 fm. 

The solid lines in Fig. | show the results of fitting data taken on a 16^ X 24 X A", lattice 
to an exponential with a constant, c + ae~^^". This fit reproduces all five data points with 
an acceptable xV^of for both M = 1.7 (xV^of = 1.1) and M = 2.1 (xV^of = 2.2). In 
contrast an exponential fit without a constant, ae~^^% does not work, resulting in a large 
xVdof > 30 - 60 (dotted lines). 

In order to examine the stability of the fit, we carry out fits to the four points excluding 
the data for the smallest fifth dimensional length A^s = 10. As seen in Table |T|, the data 
are again fitted well with an exponential plus a constant, but not without a constant. The 
constant c agrees between the five- and four-point fits within statistical errors. 
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Our data at (3 = 6.0 are shown in Fig. |^ where, for comparison, the result at (3 = 5.65 
and M = 1.7 are recapitulated from Fig. |[ One generally expects that the domain- wall 
formulation works better at weaker couplings [Q. Indeed, the decrease of m^g as a function 
of Ng is much more rapid at f3 = 6.0, giving m^g an order of magnitude smaller at /3 = 6.0 
than at P = 5.65: 

_ r 10.32(46) MeV M = 1.7 at /? = 5.65 , . 

""'^ ~ I 0.79(16) MeV M = 1.8 at /? = 6.0. ^ 

Nonetheless, we still observe a clear flattening of m^g toward a large Ng. A plain exponential 
fit without a residual constant does not work. This is true even if we drop the Ng = 4 data. 
The fit including a constant is acceptable as shown by the solid lines in Fig. |^ (see Table |I| 
for details of analyses). 

As a further test for the presence of a non-zero constant, we extract the decay rate ^ in 
Ns not only at the chiral limit {rrif = 0) but also at a finite ruf. According to the transfer 
matrix description of DWQCD [^|TU|,|TT|, the mass term, which appears only at the boundary 
in the fifth dimension, does not affect the large A^^ dependence of m^g] hence we expect that 
^ does not depend on rrif. The decay rates ^ from a fit with an exponential and a constant 
are plotted in Fig. |^, which show that this expectation is well satisfied at a strong coupling, 
and within estimated errors also at a weaker coupling. The decay rate extrapolated to the 
chiral limit from finite is consistent with the value directly extracted for rrif = 0, showing 
the consistency of the analysis. 

We conclude that a residual constant remains in the anomalous quark mass in the A^^ — >■ 
oo limit. The presence of a residual constant, even if it is small, means that quenched 
DWQCD does not realize the expected chiral symmetry for the plaquette action, at least, 
up to a-^ ~ 2 GeV. 



We remark that our conclusion differs from Ref. |]20[, which indicates that m^g vanishes 
exponentially in Ng at (3 = 6.0 with the plaquette action. We find that the two data for m^g 
are mutually consistent in the overlapping range of A^, (see Fig. ^, where open points are 
from Ref. [|20[ and the solid points denote our data). The choice of A^^ < 10 in Ref. ^| 
too small to observe the asymptotic flattening. 



was 



V. RENORMALIZATION GROUP IMPROVED ACTION 

We may ascribe the failure in fulfilling chiral symmetry seen in the last section to the 
roughness of gauge configurations with the plaquette action at the lattice spacings we have 
studied. We suspect that the massless mode may exist on sufficiently smooth gauge configu- 



rations close to the continuum limit |3i,35]. One way to realize smooth gauge configurations, 
yet keeping the lattice spacing coarse enough to make a computation feasible, is to employ 
improved gauge actions. So we discuss the case with the RG-improved gauge action in 
this section. 

The anomalous quark mass in the nominal chiral limit m^q{mf = 0,Ns,M) is plotted 
as a function of A^s in Figs. ^ and ^j. Fig. ^ compares a typical result at a strong coupling 
(/5 = 2.2) with that at a weaker coupling (/? = 2.6). Fig. ^ displays the results at P = 2.2 
for three different choices of M and for two different choices of the spatial lattice size A^-. A 
comparison of solid points (16^ x 24 x A^, lattices) and open points (12^ x 24 x A^^) verifies 



7 



the absence of a significant finite size effect at this (3. In order to take the Ng oo hmit we 
carry out exponential fits, with or without a constant, as we did for the plaquette action. 
The fit parameters are summarized in Table PTI] . 

We first discuss the results in the weak coupling region displayed in Fig. ^ As one can 
see from large values of x^/dof in Table neither c + ae~'^^* (x^/dof = 12.8) nor ae"'^^" 
(x^/dof = 16.9) fits the data well in so far as the point at A^s = 4 is included. This means 
that the leading exponential term does not dominate the two point correlators owing to a 
non-trivial overlap of quark wave functions between the two domain walls for this small Ng. 

If the data at A^^ = 4 are excluded, however, the four points for A^^ > 10 are fitted equally 
well with either c + ae'^^" (solid line) or ae~'^^° (dotted line) as seen in the figure. In fact, 
the fit with a constant yields c = —0.71(63) x 10~^, consistent with zero at one standard 
deviation; the four data points fall on ae"'^^" , indicating the realization of chiral symmetry. 

We obtain a different result for a strong coupling. Solid lines in Fig. H show fits to all five 
data points with an exponential plus a constant. They provide reasonable fits to all data. 
On the other hand, a plain exponential does not fit the data. 

Repeating the analysis for the four data points excluding A^^ = 10, we still find the 
fit with a non-zero constant (x^/dof = 0.1) being much better than a plain exponential 
(x^/dof = 4.0) at M = 2.1. The constant c is consistent between the four- and five-point 
fits. This analyses supports the conclusion that chiral symmetry is not realized at M = 2.1. 
The four data points at M = 1.7, on the other hand, are fitted well with either of the 
two forms with x^/dof < 1. The constant c for the four-point fit is consistent with zero 
within two standard deviations. Thus the possibility that m^q vanishes exponentially in A^^ 
at M = 1.7 cannot be excluded from this analysis at m/ = alone. 

To further explore this issue, we attempt to fit the data for rrif ^ 0. We find that the 
plain exponential does not fit the data at non-zero ruf, leading to a large x^/dof as shown in 
Fig. ^ Since we expect the anomalous quark mass to depend little on m/ from the transfer 
matrix formalism of DWQCD, this suggests that the good fit we obtained without a constant 
at m/ = is perhaps accidental. While further data are needed for a definitive conclusion, 
we think it likely that a non-zero anomalous quark mass also remains in the Ng ^ oo limit 
at M = 1.7. 

The value of the residual quark mass is small. From the five-point fit, we obtain 

msg = 3.60(45) MeV at M = 1.7 (V.l) 
= 2.87(23) MeV at M = 2.1 (V.2) 

We finally show the decay rates ^ at finite ruf and = in Fig. |T^. The decay rates 
at /3 = 2.2 are obtained by the five-point fit with a constant, while those at /3 = 2.6 employ 
plain exponential fits to the four data points. As expected, the decay rates depend little on 
m/; their values at m/ 7^ are consistent with those at m/ = 0. 

We conclude that the RG-improved action significantly improves the chiral behavior of 
DWQCD. While chiral symmetry is still not realized at ~ 1 GeV, the residual anomalous 
quark mass is sizably reduced compared with that with the plaquette action. The results at 
a weaker coupling (a~^ ~ 2 GeV) are consistent with the vanishing anomalous quark mass, 
supporting the realization of chiral symmetry with DWQCD. 
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VI. DISCUSSION 



Let us now try to understand our findings. The anomalous quark mass m^q vanishes 
exponentially in A^^ if the eigenvalues of the transfer matrix in the fifth direction are strictly 
less than unity []5|, p!0| , |3^ , |35| . The occurrence of a unit eigenvalue, in turn, is in one-to-one cor- 
respondence with that of a zero eigenvalue of the four dimensional Hermitian Wilson-Dirac 
operator Hy/ |^,^ or more complicated one, Hw, defined through the fifth dimensional 
transfer matrix T P,p!0|, 

It has been argued that zero eigenvalues of Hy/ are related to the existence of the parity- 



flavor broken phases in ordinary lattice QCD with the Wilson quark action The 
connection follows from the identity for the parity-flavor order parameter (qi'y^T^q) given by 



lim (qi'y^T^q 



lim Tr- 

H^+o Dw + i'^-,T^H 

1 1 

-I lim tr 



lim tr 



^75 



?75 



Dw + ilbH Dw - ilbH 



Hw + iH Hw - iH 



X + iH X~iH 
I / dXpH^{X){-2'Ki)5{X) = -27rpHH.(0), 



(VI.l) 



where H is an external field coupled to the order parameter, Hw = l^Dw with Dw the 
Wilson-Dirac operator, and PhwW is the density of the eigenvalues of Hw- 

In Fig.0, the expected phase diagram is schematically drawn in the {(3, M) plane, where 
—M is the bare quark mass of the Wilson quark action. In the shaded regions the parity- 
flavor symmetry is broken spontaneously corresponding to a non-zero density of zero eigen- 
values of Hw- At the critical lines that form the phase boundaries, the neutral pion mass 
vanishes. At a strong coupling {(3 < jSc) there are only two critical lines, while there exist 
ten of them at a weak coupling [jS > jSc) and 5 points where two lines meet at /5 = oo. Each 
point corresponds to one continuum limit, whose low energy spectra are composed of a part 
of sixteen fermion doublers. 

For DWQCD to work we have to tune the domain-wall height M within the thick shaded 
region in Fig. |ll] in order to avoid zero eigenvalues of Hw- In the weak coupling limit, this 
region is given by < A-f < 2. As the coupling increases, the range is shifted to a larger 
value of M and shrinks in width. Finally no massless fermion exists in the strong coupling 
region at P < Pc- 

A naive interpretation of our results for DWQCD according to this picture would be that 
Pc > 6.0 for the plaquette action while 2.6 > Pc > 2.2 for the RG-improved action, or in 
terms of the lattice spacing evaluated at P = Pc, ci~^>2 GeV for the former action and 
2 GeV>a-i>l GeV for the latter. 

This interpretation seems inconsistent with the quenched spectrum result for the pla- 
quette action obtained in Ref. I^Qj; in this work, non-zero pion mass, signalling the parity 
conserving phase, was reported over the range 0.909(11) < M < 2.47(5) at P = 6.0, in- 
dicating Pc < 6.0. A possible explanation reconciling this apparent conflict is as follows. 
If the order parameter, eq. ([VI. 1|) , remains non-zero in the H +0 limit, the charged 
pions become massless Nambu-Goldstone modes associated with spontaneous parity-flavor 
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breaking in the infinite volume limit. The fact that m^q in the proper chiral limit is tiny at 
P = 6.0 suggests that the magnitude of the order parameter is extremely small. Since the 
simulation in Ref. is performed on a finite volume (16^) without adding the external 



field H, it may be difficult to detect the existence of massless pions corresponding to such a 
small order parameter. The anomalous quark mass, ms^, on the other hand, is very sensitive 
to the small eigenvalues of Hw, and hence to the non-zero order parameter. To clarify this 
issue, a more detailed study for the parity-flavor breaking phase of ordinary lattice QCD 
with the Wilson quark will be needed. 

The actual dynamics inducing non-zero PHwi^) t)e complicated, possibly involving 
instantons. The relation between the spectral gap of the Hermitian Wilson-Dirac operator 



H\Y and the instanton number has been extensively studied While a single instanton 



causes a zero eigenvalue of Hw only for a single value of M, an ensemble of instantons may 
lead to a non-zero density PHwi^) 0- particular unphysical short- distance topological 
dislocations, having an action less than that of physical instantons that are present for the 
plaquette action, may lead to such an effect They will lead to an anomalous quark 

mass even at ~ 2 GeV as we have found in our simulations. 

This mechanism may also explain the difference in the lattice spacing needed to realize 
chiral symmetry in DWQCD between the plaquette and RG-improved actions. The RG- 
improved gauge action suppresses unphysical dislocations by pushing the action above that 
of physical instantons . Hence their influence on the density of zero-eigenvalues will be 
negligibly small compared to that for the plaquette action. Work exploring zero eigenvalues 
of the Wilson-Dirac operator for both gauge actions should shed light on this interesting 
problem. 



VII. CONCLUSIONS 

We have investigated the chiral property of domain-wall QCD within the quenched ap- 
proximation. We have defined the anomalous quark mass (axial WT identity quark mass) as 
an indicator for the realization of chiral symmetry. Our simulations have been made in both 
strong and weak coupling regions corresponding to ~ 1 and 2 GeV, using the plaquette 
and RG-improved gauge actions. We have found that the anomalous quark mass remains 
non-zero in the DWQCD chiral limit, — > and Ng —>■ oo for the plaquette action for 
a^^ < 2 GeV. The magnitude of chiral symmetry breaking rapidly decreases with lattice 
spacing, but the exact chiral symmetry is not realized at least for this lattice spacing. 

On the contrary, our analysis for the RG-improved action reveals a much improved 
chiral behavior. The anomalous quark mass vanishes exponentially with this action at 
~ 2 GeV, indicating that the exact chiral symmetry is realized. At ~ 1 GeV the 
improvement of the chiral behavior is not sufficient to remove a non-zero anomalous quark 
mass. 

Overall, quenched domain-wall QCD at ~ 1 GeV appears no better than usual lattice 
QCD with the ordinary Wilson quark action. Since the effect of explicit chiral symmetry 
breaking is non- negligible, the operators relevant for the electroweak matrix elements mix 
nontrivially between different chiralities. There is no a priori reason to expect that these 
mixing coefficients are small. 
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The situation is better in the weak couphng region, e.g., at ~ 2 GeV. Even with 
the plaquette action, the chiral symmetry breaking effect is significantly smaller than the 
physical u, d quark masses, and hence may not seriously affect the chiral properties of weak 
matrix elements [|l^. Moreover, the chiral property of DWQCD is distinctly improved with 
the RG action: not only the explicit breaking is shown to vanish exponentially in Ns-, but 
also the size of the breaking effect at finite A^^ is smaller. The improved gauge actions, such 
as the RG-improved one, should be a preferred choice for future numerical simulations for 
DWQCD. 
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TABLE L Number of configurations for the plaquette action and the RG improved action(RG) 
for each Ns- The anomalous quark mass is measured only at the points where number of configu- 
rations is written in boldface. 





M 


fitting range 




c + ae'^^" 










otNs 


c 


a e 




a 




xVdof 




1.7 


10-50 


0.01032(46) 


0.115(5) 0.0947(44) 


1.07 


0.087(9) 


0.0488(52) 


59.1 


5.65 


20-50 


0.00918(107) 


0.08(2) 0.077(12) 


0.00193 


0.051(7) 


0.0337(41) 


9.05 


2.1 


10-50 


0.00725 (63) 


0.061(3) 0.0838(69) 


2.23 


0.050(4) 


0.0438(46) 


34.8 




20-50 


0.00524 (149) 


0.040(6) 0.056(13) 


0.202 


0.034(2) 


0.0313(21) 


2.27 


6.0 


1.8 


4-30 


0.000565 (147) 


0.12 (1) 0.364 (22) 


8.70 


0.11 (2) 


0.333 (38) 


47.8 


1.8 


10-30 


0.000396 (79) 


0.023 (6) 0.195 (25) 


0.162 


0.012 (3) 0.123 (19) 


7.29 



TABLE IL Exponential fit of the anomalous quark mass m^q for the plaquette action. Strong 
coupling results at /3 = 5.65 are for the lattice size 16^ x 24, and those at a weak coupling of /? = 6.0 

are for 16^ x 32 lattice. Fits with all of the five date points are represented by the fitting range 
10-50 and 4-30. Four-points fits without Ng = 10 data are represented by 20-50 and 10-30. 
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p M 



fitting range 
of iV, 



c + ae~^^' 
a 



ae 



xVdof 



a 



xVdof 



1.7 



10-50 0.00371(46) 0.071(9) 0.121(12) 3.03 0.045(9) 0.066(11) 35.7 

20-50 0.00230(107) 0.025(7) 0.0620(201) 0.00002 0.020(2) 0.0370(28) 0.936 

10-50 0.00296 (24) 0.036(3) 0.0979(74) 0.232 0.025(4) 0.0505(76) 21.0 

20-50 0.00278 (43) 0.03(1) 0.085(21) 0.113 0.014(2) 0.0328(51) 3.98 



2.2 



2.1 



2.6 



1.8 4-24 0.0000099 (71) 0.122(9) 0.501(15) 
1.8 10-24 -0.0000071(63) 0.014(3) 0.281(25) 



12.8 0.12(1) 0.498(17) 16.9 
0.00227 0.019(2) 0.312 (12) 0.799 



TABLE III. Exponential fit of the anomalous quark mass m^g for the RG-improved action. 
Strong coupling results at /3 = 2.2 are for the lattice size 16^ x 24, and those at a weak coupling of 
P = 2.6 are for 16^ x 32 lattice. Fits with all of the five date points are represented by the fitting 
range 10-50 and 4-30. Four-points fit without Ng = 10 data are represented by 20-50 and 10-24. 
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APPENDIX A: PION MASS 



The pion mass is a standard quantity in examining the chiral property, and hence often 
studied in the context of DWQCD. Analyses of this observable, however, is not straightfor- 
ward due to possible finite size effects, particularly for small pion masses, and difficulties 
associated with the chiral extrapolation because of chiral logarithms, which may be quite 
significant in quenched QCD. Numerically the pion mass has to be extracted from an ex- 
ponential falloff of the pion propagator. Even though the pion propagator has the best 
statistical quality among hadron propagators, this procedure is more susceptible to statisti- 
cal and systematic uncertainties than the determination of the anomalous quark mass which 
involves only a constant fit to a ratio of two kinds of propagators. For these reasons we 
have employed the anomalous quark mass in the main body of the present paper. In this 
appendix we present our results for the pion mass, and discuss to what extent they match 
with those obtained with the anomalous quark mass. 

Typical results for pion mass squared ■m'^{mf, Ng, M) are plotted as a function of the 
averaged bare quark mass rrij" = (wiji + mj2)/2 for the plaquette action in Fig. |12[ The 
result shows a good linear behavior in mf. The pion mass in the chiral limit ruf — is 

mass are collected in Appendix B. 



therefore estimated by linearly extrapolating in mf- The numerical values of the pion 



Strong coupling region 

We first discuss pion mass results at the strong coupling a'^ ~ 1 GeV. In Fig. |l^ we plot 
the pion mass squared in the chiral limit mf = as a function of Ng for the plaquette 
action at /3 = 5.65. Filled symbols are obtained on a lattice of spatial size N^^ = 16 and 
open ones with a spatial size of N^r = 12. A similar figure for the RG-improved action is 
shown in Fig. |l^. 

Solid lines in these figures show results of fits by a form c + ae~'^^° using all five points 
obtained on an = 16 lattice (M = 1.7,2.1) or an = 12 lattice (M = 1.3,2.5). They 
reproduce the data well for a non-zero c with an acceptable ^ summarized in Tables ^ 
and 0. Therefore pion mass does not vanish in the chiral limit m/ — > and iV^ — >^ oo at 
a^^ ~ 1 GeV for the plaquette and the RG-improved actions, which is consistent with the 
conclusions from the analysis of the anomalous quark mass. 

In order to check our procedure of taking the chiral limit we interchange the order of the 



limits mY and Ng ^ oo. As shown by solid lines going through open circles in Fig. |T^, 
we first make a fit of form m^{mj' , Ng) = d{mj') + ae'^^" for each value of m^^. We note 
that the constant term has to be included in the fit of the pion mass at non-zero mj. A linear 
chiral extrapolation m'^{mJ',Ng = oo) = c'{mj") = d + 'ymj" then yields d = 0.0432(52), 
which agrees well with the value of c = 0.0440(43) previously obtained. The commutativity 
of the two limits is summarized by two symbols "X" and "Y" in Fig. |15|, where X represents 
m^{mj' = 0,Ng = oo) given with mj — > 0, then Ng ^ oo and Y represents the pion mass 
with the limit Ng oo, then mj — > 0. These two values agree very well with each other. 
The same consistency check is also made for the RG-improved action. 

The above analyses strongly support the conclusion that ttItt does not vanish in the limit 
mY and A^,, ^ oo. In order to examine finite spatial volume effects as a possible 
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origin of this non-zero mass, we compare the pion mass squared in the proper chiral hmit 
m^(m™ 0, Ng ^ oo) for two kinds of lattice volumes 12^ x 24 (open symbols) and 16^ x 24 
(filled symbols) in Fig. |16|. 

The circles and diamonds represent the results from the plaquette and the RG action, 
respectively. The horizontal triangles are results for the Nambu-Goldstone pion mass of the 
Kogut-Susskind quark action at /? = 5.7 for a spatial lattice size of N^j = 12 (open triangle) 
and Ncr = 16 (filled triangle) [^. From comparison of these data points, we consider that 
there may be some finite size effects at M = 1.7, but that they are not large enough to 
explain non-zero values of pion mass for DWQCD at strong coupling. 

Since there still remains a possibility that this non-zero mass is caused by the quenching 
effect, we adopt the WT identity mass as a better indicator to test the chiral property of 
DWQCD in the main body of this article. 



Weak coupling region 

Pion mass squared in the chiral limit m^(mj = 0, A'^5,M) is plotted in Fig. O for the 



plaquette action at /3 = 6.0 and in Fig. |18| for the RG-improved action at /? = 2.6. The 
results have large errors, and do not show a clear trend as a function of Ns. However, the 
magnitude beyond Ng = 10 are comparable to the value obtained for the Nambu-Goldstone 
pion mass of the Kogut-Susskind fermion at /3 = 6.0 for the same spatial lattice size 16^ x 40 



(open circle) P^^Sj . This indicates that these non-zero pion masses are mainly caused by 
finite spatial volume effects. Hence the chiral property of DWQCD cannot be studied by 
pion mass in the weak coupling region on this volume. 



action 


M 


fitting 


range 




c -|- ae' 












oiNs 


c 


a 


e 


xVdof 


a 




xVdof 




1.7 
2.1 


10^ 


-50 


0.04404(434) 


0.51(3) 


0.0836(62) 


0.838 


0.41(4) 


0.0485(54) 


16.1 


plaq. 


20^ 
10- 


-50 
-50 


0.03946(829) 
0.07796(464) 


0.4(1) 
0.35(3) 


0.071(16) 
0.0832(78) 


0.944 
0.761 


0.26(4) 
0.29(3) 


0.0341(46) 
0.0293(54) 


3.63 
20.4 




20- 


-50 


0.0602(262) 


0.21(5) 


0.046(28) 


0.00354 


0.21(1) 


0.0195(16) 


0.523 




1.7 


10- 


-50 


0.02549(321) 


0.30(4) 


0.109(14) 


0.585 


0.20(4) 


0.0513(92) 


8.65 


RG 


20- 


-50 


0.0191(105) 


0.14(8) 


0.063(38) 


0.0661 


0.10(1) 


0.0300(46) 


0.456 


2.1 


10- 


-50 


0.04453(426) 


0.20(4) 


0.100(21) 


0.553 


0.14(2) 


0.0285(66) 


7.17 




20- 


-50 


0.04669(449) 


0.5(1.1) 


0.15(11) 


0.677 


0.10(2) 


0.0167(59) 


2.80 



TABLE IV. Exponential fit of pion mass squared in strong coupling region at M = 1.7 and 
2.1. The four-dimensional lattice size is 16^ x 24. Fits with all of five data points are represented 
by the fitting range 10 — 50 and four-points fits without A^^ = 10 data are represented by 20 — 50. 
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action M 



fitting range 
of iV. 



a ^ 



ae 



X^/dof a 



xVdof 



plaq. 



1.3 
1.7 
2.1 

2.5 



10-50 
20-50 
10-50 

20-50 
10-50 
20-50 
10-50 

20 5(J 



0.09328(709) 
0.0842(126) 
0.0593(121) 
-0.021(281) 
0.08112(660) 
0.0765(140) 
0.14777(566) 



1.00(5) 
0.7(2) 
0.6(1) 
0.2(2) 
0.43(6) 
0.3(2) 
0.6(2) 



0.0904(57) 
0.075(14) 
0.100(25) 
0.021(56) 
0.098(14) 
0.077(42) 
0.159(33) 



0.702 
0.147 
5.27 
4.90 
0.226 
0.192 
0.617 



0.8(1) 
0.45(7 
0.40(6 
0.22(3 
0.31(4 
0.20(2 
0.29(6 



0.0500(64) 
0.0320(50) 
0.0454(74) 
0.0267(49) 
0.0332(65) 
0.0189(35) 
0.0205(86) 



27.0 
5.19 
16.6 

2.47 
11.7 
1.46 
30.2 



1.3 

RG 1.7 
2.1 



10-50 
20-50 
10-50 
20-50 
10-50 
20-50 



0.03482(494) 0.83(7) 0.1125(84) 0.437 0.63(9 



0.03175(973) 
0.03681(390) 
0.0327(152) 
0.0344(181) 
-0.035(239) 



0.6(4) 
0.4(1) 
0.1(2) 
0.13(2) 
0.2(2) 



0.098(34) 
0.150(31) 
0.08(12) 
0.051(25) 
0.013(35) 



0.670 0.31(6 

0.704 0.20(4 

1.30 0.08(2 

0.281 0.15(1 

0.1593 0.12(2 



0.0757(98) 
0.0475(69) 
0.048(11) 
0.0189(60) 
0.0262(35) 
0.0209(56) 



11.4 
2.23 
9.15 
0.759 
0.529 
0.0883 



TABLE V. Exponential fit of the pion mass squared in strong coupling region at 
M = 1.3,1.7,2.1,2.5. The four dimensional lattice size is 12^ X 24. Fits with ah of five data 
points are represented by the fitting range 10 — 50 and four-points fit without Ng = 10 data are 
represented by 20 — 50. 
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APPENDIX B 

In this appendix we lists results for anomalous quark mass m^g and pion mass squared 



nif 




10 20 30 40 


50 




P = 5.65, M = 1.3 on 12^ x 24 


0.00 




- - - 0.0605 (31) 


0.0499 (20) 


0.03 




- - - 0.0622 (28) 


0.0520 (19) 


0.05 




- - - 0.0629 (26) 


0.0530 (19) 


0.1 




- - - 0.0658 (24) 


0.0565 (21) 




13 = 5.65, M = 1.7 on 12^ x 24 


0.00 




- - - 0.01303 (56) 


0.01019 (47) 


0.03 




- - - 0.01304 (53) 


0.01019 (42) 


0.05 




- - - 0.01302 (59) 


0.01016 (42) 


0.1 




- - - 0.01302 (63) 


0.01018 (44) 




13 = 5.65, M = 2.1 on 12^ x 24 


0.00 




- - - 0.00934 (28) 


0.00770 (42) 


0.03 




- - - 0.01031 (31) 


0.00894 (29) 


0.05 




- - - 0.01051 (35) 


0.00907 (33) 


0.1 




- - - 0.01230 (55) 


0.01155 (64) 


TABLE VI. Anomalous quark mass for the plaquette action in the strong 


couphng region at 


P = 


5.65, M = 


1.3, 1.7,2.1 on 12^ x 24 lattice. 
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40 


50 






B — 


5.65, M = 1.3 on 16^ x 24 




n nn 


n 1 747 ( 1 (\) 


n 1129 ('14') 


— 


— 


n n4q25 (qr) 




n 1 72S (^ (\) 


n 1 1 23 (^Pi) 


— 


— 


n n5i '\ (^^) 


0.05 


n 1700 (17) 


1111 (n) 


— 


— 


0522 fll) 


0.1 


1660 (18) 


1102 ('20) 






0557 fl4) 






:j = 


5 (iF) Af — 1 7 on 


IG'' X 2 1 




n nn 


n n^^i4 (dQ) 


n n27i 7 ('fi9') 


0.01748 (55) 


0.01303 (27) 


n ni nQ^ (d7) 


n ns 


n n5S7s Mfi^i 


n n2646 (f>R) 


0.01721 (51) 


0.01280 (29) 


n ni n92 c^q) 


0.05 


0.05228 (49) 


0.02561 (59) 


0.01691 (51) 


0.01252 (35) 


0.01084 (43) 


0.1 


0.05002 (67) 


0.02445 (60) 


0.01645 (54) 


0.01219 (52) 


0.01080 (44) 






/? = 


5.65, M = 2.1 on 


16^ X 24 




0.00 


0.03349 (31) 


0.01820 (37) 


0.01253 (39) 


0.00954 (24) 


0.00757 (43) 


0.03 


0.03338 (27) 


0.01941 (35) 


0.01358 (37) 


0.01032 (22) 


0.00880 (38) 


0.05 


0.03261 (28) 


0.01947 (38) 


0.01361 (37) 


0.01032 (24) 


0.00909 (46) 


0.1 


0.03254 (37) 


0.02171 (56) 


0.01542 (44) 


0.01203 (56) 


0.01124 (51) 



TABLE VII. Anomalous quark mass for the plaquette action in the strong couphng region at 
(3 = 5.65, M = 1.3, 1.7,2.1 on 16^ x 24 lattice. 
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0.00084 (11) 


0.000467 (57) 


0.02 


0.02884 (25) 


0.00360 (11) 


0.00134 (14) 


0.000750 (90) 


0.000403 (49) 


0.04 


0.02804 (20) 


0.003446 (78) 


0.001128 (93) 


0.000637 (73) 


0.000317 (42) 


0.06 


0.02735 (18) 


0.003316 (64) 


0.000989 (66) 


0.000548 (59) 


0.000259 (38) 



TABLE VIII. Anomalous quark mass for the plaquette action in the weak coupling region at 
p = 6.0, M = 1.8 on 16^ x 32 lattice. 



20 



7T7 f 
IILf 




10 20 30 40 


50 




/? = 2.2, M = 1.3 on 12^ x 24 


0.00 

0.03 
0.05 
0.1 




- - - 0.01565 (90) 

- - - 0.01605 (84) 

- - - 0.01657 (86) 

- - - 0.0172 (11) 


0.01260 (81) 
0.01324 (78) 
0.01385 (80) 
0.01491 (84) 




3 = 2.2, AI = 1.7 on 12'^ x 24 


0.00 
0.03 
0.05 
0.1 




- - - 0.00413 (23) 

- - - 0.00399 (20) 

- - - 0.00403 (25) 

- - - 0.00377 (21) 


0.00319 (26) 
0.00302 (24) 
0.00286 (25) 
0.00257 (26) 




p = 2.2, M = 2.1 on 12^ x 24 


0.00 

0.03 
0.05 
0.1 




- - - 0.00378 (25) 

- - - 0.00386 (21) 

- - - 0.00373 (22) 

- - - 0.00388 (23) 


0.00303 (24) 
0.00306 (19) 
0.00299 (17) 
0.00304 (24) 


TABLE IX. Anomalous quark mass for the RG improved action in the strong 
at ^ = 2.2, M = 1.3, 1.7, 2.1 on 12^ x 24 lattice. 


coupling region 



21 





10 


20 


Ns 

30 


40 


50 






B — 


2.2, M = 1.3 on 16^ x 24 




n nn 


n nQ72Q (QQ) 




0.02340 (41) 




m 329 (4.R) 




n nQ,^ (^ ^ 




0.02345 (39) 




n m 356 ('49') 


0.05 


0928 ilV] 




0.02340 (41) 




01378 (59) 


0.1 


0890 fill 




0.02350 (55) 




01422 f75) 






:-] = 


2.2, j\I = 1.7 on 1 


X 24 




n nn 


n n948fi (A')\ 


n nnQfifi fs^^i 


0.00626 (39) 


0.00443 (28) 


n nns44 fss"! 


n ns 


n n2S87 

\J . L/^OO 1 I (JO J 


n nnq2n c^?) 


0.00586 (35) 


0.00412 (26) 


n nn329 (9R) 


0.05 


0.02281 (35) 


0.00875 (33) 


0.00549 (33) 


0.00379 (23) 


0.00321 (22) 


0.1 


0.02113 (35) 


0.00800 (32) 


0.00480 (33) 


0.00326 (24) 


0.00296 (20) 






P = 


2.2, M = 2.1 on 16^ x 24 




0.00 


0.01633 (35) 


0.00788 (30) 


0.00490 (26) 


0.00377 (25) 


0.00314 (22) 


0.03 


0.01626 (29) 


0.00804 (28) 


0.00511 (21) 


0.00381 (19) 


0.00318 (18) 


0.05 


0.0158 (3) 


0.00782 (29) 


0.00501 (20) 


0.00368 (18) 


0.00312 (16) 


0.1 


0.01570 (28) 


0.00811 (30) 


0.00535 (21) 


0.00374 (20) 


0.00317 (19) 



TABLE X. Anomalous quark mass for the RG improved action in the strong couphng region 
at ^ = 2.2, M = 1.3, 1.7, 2.1 on 16^ x 24 lattice. 



ruf 






Ns 








i 


10 


16 


20 


21 


0.0 


0.01649 (11) 


0.000817 (20) 


0.000146 (22) 


0.000042 (11) 


0.000009 (2) 


0.02 


0.016238 (99) 


0.000812 (17) 


0.000136 (17) 


0.0000405 (87) 


0.0000087 (15) 


0.04 


0.015920 (80) 


0.000803 (16) 


0.000116 (12) 


0.0000397 (78) 


0.0000090 (19) 


0.06 


0.015677 (70) 


0.000801 (16) 


0.0001060 (76) 


0.0000387 (73) 


0.0000089 (19) 



TABLE XL Anomalous quark mass for the RG improved action in the weak coupling region 
at /? = 2.6, M = 1.8 on 16^ x 32 lattice. 



22 



m/ 


1 A 




20 




Ns 

30 




40 




50 












13 = 


5.65, M = 1.3 on 12^ x 24 








0.00 


A /I A O O 

0.4982 


/' oa\ 

(80) 


0.2504 


(83) 


0.1618 


(47) 


0.1234 


(71) 


0.1007 


(59) 


U.U6 


0.5ol7 


(81) 


A O 1 A'7 

0.3197 


/ o o\ 

(82) 


0.2258 


(45) 


0.1845 


(66) 


0.1597 


(61) 


n n A 


O.o07o 


(81) 


0.3415 


/ o o\ 

(82) 


0.2462 


(44) 


0.2038 


0\ 

(63) 


A 1 TOT 

U.iYOY 


id A \ 

(64) 


U.Uo 


0.D341 


/ oo\ 

(82) 


0.3635 


(83) 


0.2666 


(43) 


0.2233 


I'd \ 
(61) 


A 1 ATO 


(66) 


U.Uoo 


0.d7o0 


(83) 


0.3981 


f O A\ 

(84) 


0.2987 


(44) 


0.2536 


(61) 


A OOTO 


(71) 


[).[)( 


A 'VAO A 

0.7U3U 


(83) 


A A OA£? 

0.4206 


(84) 


0.3195 


(43) 


0.2736 


(59) 


A O /I C/^ 

(J.z4bD 


(72) 


O.i 


0.7733 


(84) 


A A 'ynci 

0.4790 


(85) 


0.3733 


(43) 


0.3249 


(59) 


A OA£? A 

0.2964 


(76) 










P = 


5.65, M = 1. 


7 on 


12^ X 24 








A AA 

0.00 


0.2d75 


(55) 


A 1 OAO 

0.1292 


(68) 


0.0889 


(61) 


0.0833 


(55) 




(50) 


n no 
O.Oo 


0.4190 


(55) 


0.2840 


(61) 


0.2440 


(53) 


0.2384 


(54) 


A 01 1 A 

U.zll4 


(50) 


A A /I 

U.U4 


0.4dd7 


(57) 


A OOOA 

0.3329 


(61) 


0.2927 


(53) 


0.2882 


(54) 


A O/^ O A 

0.2630 


(54) 


A AC 

U.Uo 


A 1 /I O 

0.5148 


(59) 


A O OOO 


(61) 


0.3416 


(52) 


0.3375 


(52) 


A O 1 /I 7 

0.3147 


(58) 


A A/^ 


A r" A 1 o 

0.5912 


(62) 


A /I CT OA 

0.4589 


(63) 


0.4193 


(55) 


0.4171 


(61) 


A OAO A 

0.3929 


(64) 


0.075 


A /I A1 

O.o401 


(d3) 


A CAA A 

0.5094 




0.4692 


(54) 


0.4662 


A\ 

(59) 


r\ A A CO 

0.4453 


(65) 


A 1 

0.1 


0.7d79 


(o5) 


A /I A1 

0.6401 


(61) 


0.6004 


(56) 


0.5966 


(61) 


0.5781 


(69) 










P = 


5.65, M = 2. 


1 on 


12^ X 24 








0.00 


0.2418 


(68) 


0.1393 


(61) 


0.1067 


(60) 


0.0876 


(68) 


A A O O A 

0.0839 


(67) 


A AO 

O.Uo 


A A '\ A'^ 

0.4143 


(70) 


A OO A/^ 

0.3206 


(56) 


0.2903 


(55) 


0.2719 


/ r* A\ 

(64) 


A O/^ O A 

0.2689 


(54) 


A A /I 

0.04 


A A TOO 

U.47zz 


(74) 


A O O 1 O 

0.3813 


(57) 


0.3521 


(55) 


0.3344 


(65) 


A O OOT 

0.3327 


(67) 


0.05 


n cone 
(J. 5^95 


(77) 


n A A'\ Cl 

0.4416 


(57) 


0.4134 


(55) 


0.3962 


(64) 


n once c 
0.3955 


/£;a\ 

(oyj 


0.0d5 


(j.blYO 


/ oo\ 

(82) 


0.5328 


(62) 


0.5056 


(60) 


0.4895 


(69) 


0.4888 


(77) 


0.0r5 


A Ar\ 

0.D74U 


(83) 


A C AOA 

0.5929 


/ r'r\\ 

(60) 




OO J 




(67) 


A C C A C 

0.5505 


[lb) 


0.1 


0.8172 


(85) 


0.7439 


(58) 


0.7195 


(57) 


0.7037 


(68) 


A Tf\ O T 

0.7037 


(77) 










/5 = 


5.65, M = 2. 


5 on 


12^ X 24 








0.00 


0.2624 


(63) 


0.1701 


(48) 


0.1546 


(43) 






0.1430 


(83) 


0.03 


0.4007 


(59) 


0.3049 


(49) 


0.2796 


(43) 






0.2439 


(77) 


0.04 


0.4496 


(59) 


0.3535 


(52) 


0.3270 


(46) 






0.2894 


(80) 


0.05 


0.4976 


(59) 


0.4012 


(56) 


0.3729 


(49) 






0.3320 


(83) 


0.065 


0.5682 


(61) 


0.4698 


(61) 


0.4385 


(54) 






0.3880 


(88) 


0.075 


0.6146 


(60) 


0.5153 


(63) 


0.4806 


(56) 






0.4236 


(89) 


0.1 


0.7284 


(60) 


0.6239 


(68) 


0.5783 


(60) 






0.4984 


(88) 



TABLE XII. Pion mass squared for the plaquette action in the strong couphng region at 
/? = 5.65, M = 1.3, 1.7, 2.1, 2.5 on 12^ x 24 lattice. 



23 



m/ 










Ns 










10 




20 




30 


40 


50 






13 = 5.65, M = 1.3 on 16^ x 24 


0.00 


0.499 (5) 


0.2556 


(40) 






0.0957 


(53) 


U.U6 


0.5823 


(52) 


0.3257 


(42) 






0.1545 


(50) 


n n A 


0.6084 


(53) 


0.3480 


(43) 






0.1734 


(49) 


U.Uo 


0.6347 


(54) 


0.3704 


(45) 






0.1925 


(48) 


U.Uoo 


0.6766 


(55) 


0.4056 


(47) 






0.2215 


(49) 


[).[)( 


0.7035 


(56) 


0.4284 


(48) 






0.2409 


(48) 


O.i 


0.7737 


(58) 


0.4876 


(50) 






0.2905 


(46) 




/? = 5.65, M = 1.7 on 16^ x 24 


A AA 

0.00 


0.2661 


(49) 


0.1358 


(56) 


0.0888 (36) 


0.0599 (38) 


0.0520 


(31) 


n no 
O.Oo 


0.4179 


(51) 


0.2916 


(55) 


0.2434 (34) 


0.2183 (32) 


0.2086 


(30) 


A A /I 

U.U4 


0.4656 


(53) 


0.3411 


(57) 


0.2929 (37) 


0.2693 (33) 


0.2599 


(35) 


A AC 

U.Uo 


0.5137 


(55) 


0.3911 


(59) 


0.3426 (40) 


0.3204 (33) 


0.3113 


(38) 


U.Udo 


0.5898 


(58) 


0.4687 


(65) 


0.4198 (46) 


0.3992 (37) 


0.3894 


(47) 


0.075 


0.6392 


(59) 


0.5198 


(64) 


0.4707 (48) 


0.4512 (36) 


0.4414 


(48) 


0.1 


0.7676 


(61) 


0.6513 


(66) 


0.6019 (53) 


0.5849 (36) 


0.5735 


(54) 










f3 = 5.65, M = 2.1 on 16^ 


X 24 






0.00 


0.2317 


(40) 


0.1434 


(33) 


0.1132 (72) 


0.0935 (50) 


0.0815 


(38) 


0.03 


0.4052 


(46) 


0.3247 


(37) 


0.2973 (62) 


0.2792 (45) 


0.2675 


(43) 


0.04 


0.4629 


(50) 


0.3860 


(42) 


0.3604 (63) 


0.3416 (49) 


0.3317 


(47) 


0.05 


0.5202 


(55) 


0.4464 


(46) 


0.4226 (62) 


0.4034 (51) 


0.3945 


(51) 


0.065 


0.6085 


(62) 


0.5390 


(52) 


0.5166 (68) 


0.4979 (58) 


0.4902 


(60) 


0.075 


0.6652 


(65) 


0.5982 


(56) 


0.5774 (63) 


0.5589 (56) 


0.5509 


(61) 


0.1 


0.8095 


(72) 


0.7484 


(60) 


0.7299 (60) 


0.7131 (54) 


0.7036 


(66) 



TABLE XIII. Pion mass squared for the plaquette action in the strong couphng region at 
f3 = 5.65, M = 1.3, 1.7, 2.1 on 16^ x 24 lattice. 



24 



m/ 


1 A 


OA 




Ns 

30 




40 




cn 








P = 


A A 

2.2, 


M = 1.3 on 12^ X 24 








0.00 


0.3038 (55) 


0.1208 (55) 




0.0684 


^71) 


0.0415 


/ cr a\ 
(50) 


A A O O A 

0.0389 


(60) 


U.U6 


A A A AO / 0\ 

0.4002 (52 j 


0.2106 (52) 




0.1579 


^64) 


0.1325 


(47) 


A 1 A CT O 

0.1253 


(57) 


n n A 


0.4301 (52j 


AOOOO /CO\ 

0.2388 (53) 




0.1867 


^63) 


0.1619 


(49) 


A 1 C 1 O 


(59) 


U.Uo 


0.4d04 (52j 


0.2671 (54) 




0.2154 


[61) 


0.1916 


(49) 


A -I ^OCi 


(60) 


U.Uoo 


0.50o3 (52) 


0.3120 (56) 




0.2605 


(64) 


0.2363 


(55) 




(a A \ 
(64) 


[).[)( 


0.5395 (52) 


0.3410 (56) 




0.2894 


[61) 


0.2667 


(53) 


A A /I A /I 

(J.z4y4 


(63) 


O.i 


0.0209 (51) 


0.4169 (56) 




0.3647 


[58) 


0.3433 


(54) 


A OOAO 


(60) 






P = 


2.2, M = 1.7 


on 


12^ X 24 








A AA 

0.00 


A 10A/1 ('TAX 

0.1304 (50) 


0.0575 (50) 




0.0405 


[57) 


0.0408 


(42) 


A AOAA 

0.0302 


(73) 


n no 
O.Oo 


0.2892 (51) 


0.2186 (46) 




0.2013 


[54) 


0.2049 


[6\)) 


A 1 OAO 

0.1898 


(57) 


A A /I 

U.U4 


0.3394 (54) 


A A£? A*? / A 0\ 

0.2697 (48) 




0.2529 


[56) 


0.2588 


( A A^ 
(4Uj 


A A /I 1 A 

(J.Z41U 


(55) 


A AC 

U.Uo 


A OAAA / r:'7\ 

0.3900 (57) 


0.3210 (49) 




0.3042 


[58) 


0.3119 


( A A^ 
(4Uj 


A OA1 'I 

0.2917 


(51) 


U.Udo 


0.4d95 (d4) 


A Ani /cc\ 

0.4016 (55) 




0.3850 


[66) 


0.3961 


( AO\ 


A O 'TA A 

0.3724 


(56) 


0.075 


0.5213 (d4) 


0.4540 (54) 




0.4372 


[68) 


0.4486 


( ACi.\ 

(4oj 


A A OO A 

0.4234 


(53) 


0.1 


0.b559 (b7) 


0.5904 (56) 




0.5737 


[77) 


0.5863 


(47) 


0.5582 


(58) 






P = 


2.2, 


M = 2.1 


on 


12^ X 24 








0.00 


0.1159 (64) 


A A01A //^OX 

0.0810 (63) 




0.0645 


[68) 


0.0556 


(73) 


0.0415 


(82) 


A AO 

O.Uo 


0.2879 (60) 


0.2510 (55) 




0.2303 


[64) 


0.2272 


(71) 


A A 1 O 1 

0.2131 


(73) 


A A /I 

0.04 


0.3461 (63) 


0.3083 (56) 




0.2873 


[66) 


0.2848 


(75) 


A A '7 A C 

0.2705 


(74) 


0.05 


0.4038 (65) 


A A A it^O\ 

0.3644 (58) 




0.3429 


[68) 


0.3415 


(78) 


0.3274 


(72) 


0.0d5 


0.4911 (72) 


0.4519 (64) 




0.4276 


[76) 


0.4292 


(87) 


C\ AA AO 

0.4143 


(81) 


0.075 


r\ v AO { t\ 

0.548 (7) 


0.5068 (66) 




0.4817 






(85j 


A A 'yr\f* 

0.4706 


(76) 


A 1 

0.1 


0.6906 (72) 


0.6482 (71) 




0.6198 


[69) 


0.6279 


(o4j 


0.6136 


(76) 






P = 


2.2, M = 2.5 


on 


12^ X 24 








0.00 


0.1549 (48) 


0.102 (6) 




0.1127 


[63) 










0.03 


0.2814 (47) 


0.2208 (57) 




0.2184 


[53) 










0.04 


0.3256 (48) 


0.2649 (57) 




0.2611 


[54) 










0.05 


0.3692 (51) 


0.3076 (58) 




0.3018 


[53) 










0.065 


0.4329 (56) 


0.3699 (59) 




0.3607 


[54) 










0.075 


0.4754 (59) 


0.4099 (59) 




0.3965 


[52) 










0.1 


0.5793 (65) 


0.5058 (58) 




0.4797 


[48) 











TABLE XIV. Pion mass squared for the RG improved action in the strong couphng region at 
/? = 2.2, M = 1.3, 1.7, 2.1, 2.5 on 12^ x 24 lattice. 
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40 
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/J 

p = 


O O 

2.2, 


M = 1.3 on 16^ X 24 






0.00 


0.3126 (54) 






0.0709 


^19) 


— 


0.0464 


(30) 


n no 

O.Oo 


0.4095 (56 j 






0.1609 


^19) 


— 


0.1348 


(28) 


U.U4 


0.4400 (58 j 






0.1900 


^19) 


— 


0.1634 


/ oA^ 
[2\)) 


(J. (Jo 


0.470/ (59 j 






0.2190 


[21) 


— 


0.191o 


[6{j) 


(j.Uoo 


0.5194 (6zj 






0.2642 


(25) 


— 


0.2368 


A\ 

(34) 




0.5508 (boj 






0.2938 


[26) 


— 


A O/^ C O 


O A\ 

(34) 


O.i 


O.bozb (b6j 






0.3702 


[33) 


— 


A OOAC 

0.3395 


(35) 






p = 


2.2, M = 1.7 


on 


16^ X 24 






0.00 


0.1276 (49j 


A ACTA ('0'v\ 

U.Uo79 {61} 




0.0392 


[36) 


0.0308 (35) 


A AO A C 

0.0246 


(35) 


n no 
O.Oo 


0.28O2 (46) 


A OOOO ( O A\ 

[j.2222 [64:) 




0.2012 


[32) 


0.1959 (33) 


0.1908 


A\ 

(34) 


U.U4 


A O O AO / /I C \ 


0.2763 (36) 




0.2527 


[35) 


0.2498 (37) 


A O /I /I 

(J.z454 


[6b) 


U.Uo 


A OAAFC ( A A \ 

U.oyUo (44j 


0.3300 (3o) 




0.3042 


[37) 


0.3034 (40) 


(J.zyyo 


( oo\ 
[oo) 


O.Odo 


A /I '7 A /I / /I \ 

0.4704 (46) 


0.4131 (44) 




0.3854 


[42) 


0.3874 (46) 


A 0*7 

0.3837 


( A A\ 

(44j 


0.075 


0.5228 (45) 


0.4667 (45) 




0.4381 




n 4409 ('48') 


0.4378 


( Ati\ 

(45j 


0.1 


O.b5o2 (45) 


n ctr\Aci f AO\ 
0.6046 (48) 




0.5760 


[47) 


0.5798 (53) 


0.5773 


(47) 






P = 


2.2, 


M = 2.1 


on 


16^ X 24 






0.00 


0.1194 (51) 


0.0735 (40) 




0.0539 


[46) 


0.0455 (41) 


0.0492 


(46) 


0.03 


0.2868 (42) 


0.2442 (38) 




0.2303 


[45) 


0.2177 (40) 


0.2238 


(44) 


0.04 


0.3425 (40) 


0.302 (4) 




0.2912 


[49) 


0.2754 (43) 


0.2840 


(48) 


0.05 


0.398 (4) 


0.3590 (41) 




0.3509 


[51) 


0.3324 (45) 


0.3427 


(49) 


0.065 


0.4826 (43) 


0.4459 (47) 




0.4413 


[59) 


0.4202 (53) 


0.4336 


(58) 


0.075 


0.5376 (43) 


0.5019 (48) 




0.4991 


[58) 


0.4763 (51) 


0.4899 


(56) 


0.1 


0.6770 (49) 


0.6434 (54) 




0.6445 


[59) 


0.6197 (51) 


0.6336 


(56) 



TABLE XV. Pion mass squared for the RG improved action in the strong coupHng region at 
13 = 2.2, M = 1.3, 1.7, 2.1 on 16^ x 24 lattice. 
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fUf 






N 








A 

'-X. 


10 


1 6 


90 


30 


0.0 


0.0891 (49) 


0.0096 (36) 


0.0025 (27) 


0.0063 (21) 


0.0019 (20) 


0.02 


0.1617 (43) 


0.0795 (30) 


0.0711 (24) 


0.0736 (19) 


0.0703 (17) 


0.03 


0.1959 (41) 


0.1130 (34) 


0.1032 (25) 


0.1051 (20) 


0.1029 (21) 


0.04 


0.2310 (39) 


0.1471 (32) 


0.1359 (24) 


0.1371 (20) 


0.1362 (20) 


0.04 


0.2311 (40) 


0.1471 (36) 


0.1362 (26) 


0.1380 (21) 


0.1366 (23) 


0.05 


0.2671 (38) 


0.1823 (35) 


0.1710 (27) 


0.1714 (23) 


0.1706 (20) 


0.06 


0.3041 (36) 


0.2179 (31) 


0.2059 (27) 


0.2061 (23) 


0.2059 (20) 



TABLE XVI. Pion mass squared for the plaquette action in the weak couphng region at /3 = 6.0, 
M = 1.8 on 16^ X 32 lattice. 
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Ns 








4 


10 


16 


20 


24 


0.0 


0.0536 (26) 


0.0044 (33) 


0.0052 (23) 


0.0023 (19) 


0.0022 (16) 


0.02 


0.1188 (24) 


0.0670 (30) 


0.0671 (23) 


0.0646 (15) 


0.0652 (14) 


0.03 


0.1505 (26) 


0.0967 (32) 


0.0969 (26) 


0.0942 (17) 


0.0959 (15) 


0.04 


0.1826 (29) 


0.1268 (30) 


0.1270 (27) 


0.1243 (17) 


0.1268 (15) 


0.04 


0.1835 (29) 


0.1274 (34) 


0.1280 (29) 


0.1246 (19) 


0.1274 (16) 


0.05 


0.2150 (36) 


0.1583 (31) 


0.1582 (31) 


0.1559 (18) 


0.1580 (16) 


0.06 


0.2491 (39) 


0.1907 (33) 


0.1904 (31) 


0.1877 (18) 


0.1908 (16) 



TABLE XVII. Pion mass squared for the RG improved action in the weak coupling region at 
(3 = 2.6, M = 1.8 on 16^ x 32 lattice. 
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FIGURES 
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FIG. 1. Ratio of two point functions < J5g(t)P(0) > / < P{t)P{0) > for plaquette action at 
P = 5.65, M = 1.7 on a 16^ x 24 x 30 lattice (left) and at /3 = 6.0, M = 1.8 on a 16^ x 32 x 30 
(right). Lines show constant fit over the fitted range. 



RG p=2.2 16 x24x30 M=1.7 



g 0.005 



Tiii 
II 



0.000 
0.010 



g 0.005 



0.000 
0.010 



£ 0.005 



1 



I 1 5 ^ ^ ^ 5 isfiSESf^ |j| 5 X 5 (|) ;| goQQ^ 



10 15 
t 



RG improved, p=2-6, 16 x32x24, M=1.8 



O m,=0.06 



)j _ A 0.0001 



0.0002 - ^ 




FIG. 2. Ratio of two point functions < J5q{t)P{0) > / < P{t)P{0) > for RG improved action 
at P = 2.2, M = 1.7 on a 16^ x 24 x 30 lattice (left) and at ^ = 2.6, M = 1.8 on a 16^ x 32 x 24 
lattice (right). Lines show constant fit over the fitted range. 
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chiral extrapolation of anomalous quark mass 

Plaquette p=5.65 16'x24xN, M=1.7 T-range:8-16 
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Plaquette p=6.0 16 x32xN, M=1.l 
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FIG. 3. Anomalous quark mass m^q as a function of m/ for the plaquette action at M = 1.7 
on a 16'^ X 24 X A'^s lattice in the strong coupling region (left) and at M = 1.8 on a 16^ x 32 x 
lattice in the weak coupling region (right). 
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FIG. 4. Anomalous quark mass m^q as a function of in mj — > limit for the plaquette action 
at /3 = 5.65. Filled symbols are data on 16^ x 24x A^^ lattice, and open ones on 12^ x 24 x A^g lattice. 
Lines are fits to all the filled points with two functions: ae~^^'' (dotted line) and c + ae~^^° (solid 
line). 
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plaquette 



• P=5.65 (16'x24xNJ 
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FIG. 5. Anomalous quark mass m^q as a function of A'^, in m/ — > limit for the plaquette action 
at weak coupling /3 = 6.0 with M = 1.8 (filled squares) on a 16^ x 32 x A^^ lattice. Representative 
strong coupling results taken at /3 = 5.65 and M = 1.7 are reproduced from Fig. ^ for comparison 
(filled and open circles). The results at /? = 6.0 from a previous simulation |^0| is also given 
with open diamonds. Lines are fits to filled symbols with two functions: oe"^^' (dotted line) and 
c + ae"^-^" (solid line) . 
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FIG. 6. Decay rate ^ of the anomalous quark mass from the fit with c + ae"^^" for plaquette 
action. Open circles represent the data in the strong coupling region at Af = 1.7, /? = 5.65 on 
16^ X 24 X A^s lattice with five points fit. Open squares are those in the weak coupling region at 
M = 1.8, P = 6.0 on 16^ x 32 x A^^ lattice with four points fit. 
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FIG. 7. Anomalous quaxk mass as a function of iV^ in m/ ^ limit for the RG improved 

action. Filled circles represent data in the strong coupling region at M = 1.7 and (3 = 2.2 on 
16^ X 24 X A'^s, lattice. Filled squares are those in the weak coupling region at M = 1.8 and (3 = 2.6 
on 16^ X 32 X A^s lattice. Lines are fits with two functions: ae^^^^ (dotted line) and c + ae"^-^" 
(solid line). Data for Ng = 10 — 50 (5 points) and 10 — 24 (4 points) are used for the fit in the 
strong and weak coupling. 
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FIG. 8. Anomalous quark mass m^g as a function of A^s in mj ^ limit for the RG-improved 
action at /? = 2.2. Filled symbols are data on 16^ x 24 x iVg lattice, and open ones on 12^ X 24 x 
lattice. Lines are fits to all the filled points with two functions: ae~^^' (dotted line) and c+ae~^^" 
(solid line). 
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FIG. 9. x^/d-of in the five and four points fit of m^g data in strong coupling RG action at 
f3 = 2.2, M = 1.7. Open circles represent those from the four points fit with simple exponential 
form ae~^^" . 
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FIG. 10. Decay rate ^ of the anomalous quark mass at each ruf for the RG- improved action. 
Open circles represent data in the strong coupling region from the five points fit with c + ae"^^" 
at M = 1.7. The open squares are those in the weak coupling region from the four points fit with 
ae-^^- at M = 1.8. 
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FIG. 11. Phase structure of the four dimensional Wilson fermion system. M represents the 
fermion mass (domain- wall height). The five cusps at /? = oo correspond to the points where 
subsets of 16 fermions become massless. The shaded area marked by < vr represents the 
parity broken phase. 



chiral extraporation 

Plaquette M-65 16\24xN, M=1.7 T-range:7-12 
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FIG. 12. Pion mass squared as a function of valence quark mass rrij at M = 1.7, /3 = 5.65 on 
16^ X 24 X Ns lattice (left) and f3 = 6.0 on 16^^ x 32 x A^s lattice (right) of plaquette action. Lines 
are linear fits to those data. 
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FIG. 13. Pion mass squared at m/ = as a function of the extra dimension Ng in the strong 
couphng region /3 = 5.65 for plaquette action. Filled symbols are taken on 16^ x 24 x iV^ lattice, 
and open ones on 12^ x 24 x A^^ lattice. Lines are fits to an exponential with a constant of all 
five data points on an A^^- = 16 lattice except for M = 1.3 and 2.5 for which data fitted are for 

= 12. 
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FIG. 14. Pion mass squared at mj = as a function of the extra dimension Ng in the strong 
coupling region /3 = 2.2 for RG-improved action. Filled symbols are taken on 16^ X 24 X iVs lattice, 
and open ones on 12^ x 24 x A^^ lattice. Lines are fits to a constant and an exponential of all five 
data points on an A/g. = 16 lattice except for M = 1.3 for which data fitted are for Ncr = 12. 
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vs 

Plaquette p=5.65 16''x24xN, M=1.7 T-range:7-12 
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RG |3=2.2 16'x24xN, M=1.7 T-range;7-12 
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FIG. 15. Pion mass squared as a function of the extra dimension A'^s at M = 1.7 on 16^ x 24 x A'^s 
lattice. The left figure represents the results for the plaquette action at /? = 5.65. The right figure 
is for the RG action at /3 = 2.2. The pion masses derived with two different order of limits m/ ^ 0, 
Ns ^ oo are represented with "X" and "Y". 
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FIG. 16. Pion mass squared as a function of domain- wall height M in the chiral limit m/ ^ 0, 
Ns — > oo. The results on 16^ x 24 lattice (open symbols) are compared with those on 12^ X 24 
(filled symbols). Open and filled triangles are those from the KS fermion action at N„ = 12 and 
Nu = 16 lattice volumes, which represent the pure finite volume effect. 
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FIG. 17. Pion mass squared at m/ = as a function of Ng at P = 6.0 of plaquette action with 
M = 1.8 (filled diamonds) as compared with those at /3 = 5.65 and and M = 1.7 (squares). Open 
circle represents the pion mass from the KS fermion at N^r = 16 lattice volume. 
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FIG. 18. Pion mass squared at m/ = as a function of Ng at /3 = 2.6 of RG-improved action 
with M = 1.8 (filled diamonds) as compared with those at /3 = 2.2 and and M = 1.7 (squares). 
Open circle represents the pion mass from the KS fermion at N^- = 16 lattice volume. 
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